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1.  INTRODUCTION 

—The  purpose  of  this 'SBIR^program  is  to  research  and  develop  a  new  laser  Doppler 
velocimeter  for  deployment  in  the  ocean.  The  system  must  be  programmable  and 
autonomous,  and  must  operate  unattended  over  periods  of  about  six  months.  The  system 
will  ultimately  operate  in  highly  dynamic  shallow  water  environments  and  must  be  capable 
of  providing  magnitude  and  direction  of  the  velocity  distribution  in  two  directions  and 
associated  flow  turbulence  as  a  function  of  time.  The  system  will  be  utilized  to  measure 
sediment  motion  and  a  variety  of  ocean  transport  phenomena  in  highly  dynamic  situations. 

The  choice  of  hardware  was  suggested  in  the  Phase  I  proposal  and  has  been  implemented 
during  the  Phase  I  work.  It  included  a  laser  diode  to  illuminate  the  particles,  three  Bragg 
cells  to  obtain  velocity  directionality  and  to  separate  the  horizontal  and  vertical  velocity 
components,  an  avalanche  photodiode  (APD)  to  collect  the  scattered  light,  a 
programmable  digital  signal  processor  (DSP)  and  solid  state  memory  to  record  and  process 
the  velocity  data,  and  optics  to  produce  the  measurement  volume  and  collect  the  scattered 
light.  All  of  these  components  were  carefully  chosen  with  the  goal  of  satisfying  the 
NAVY’s  requirements  with  the  smallest  possible  system. 

As  indicated  later  in  this  report,  the  total  power  consumption  of  the  final  system  would  be 
about  10  watts  during  data  acquisition,  and  0.2  watts  during  the  waiting  periods.  Total  time 
of  data  acquisition  during  the  six  month  period  is  expected  to  be  about  24  hours.  This 
operational  time  assumes  that  a  velocity  reading  is  made  every  hour  on  the  hour.  This  time 
could  be  reduced  if  instead  of  collecting  data  at  constant  time  intervals,  a  sensor  such  as  a 
microphone  or  an  accelerometer  was  used  to  activate  the  velocimeter  during  times  of 
vigorous  activity.  The  data  acquisition  logic  will  be  decided  in  close  consultation  with 
NAVY  technical  personnel. 

During  the  performance  of  the  Phase  I  work  we  compared  the  capabilities  of  a  pulsed  laser 
diode  configuration  with  those  of  a  continuous  wave  laser  diode  and  concluded  that  no 
major  advantages  were  obtained  by  pulsing  the  laser  diode.  This  conclusion  would  have 
been  different  if  a  different  family  of  lasers  had  been  considered  (such  as  YAG).  However 
the  power  limitations  imposed  by  a  battery  operated  system  clearly  define  the  laser  diode 
as  the  best  candidate  ana  we  therefore  chose  it  for  our  work. 

The  velocity  of  the  water  flow  was  measured  in  a  test  cell  oriented  both  horizontally  and 
tilted  34  degrees  from  the  horizontal.  The  vertical  and  horizontal  velocity  components 
were  simultaneously  measured  thus  providing  the  magnitude  of  the  velocity  as  well  as  the 
direction  traveled  by  each  particle  crossing  the  probe  volume.  The  light  scattered  by  these 
particles  was  collected  by  lenses  and  directed  to  an  avalanche  photodiode  where  it  was 
converted  into  an  electronic  signal.  This  electronic  signal  was  amplified  and  stored  in  solid 
state  memory  for  after-processing.  FFTs  were  performed  by  the  programmable  digital 
signal  processor  providing  the  frequency  spectrum  which  separated  the  two  components  of 
the  velocity.  The  frequency  resolution  using  a  100  MHz  A/D  converter  was  estimated  at 
2.5  KHz  which  corresponded  to  0.07  m/s.  For  a  velocity  of  7  m/s  this  represents  a 
resolution  of  1%.  The  frequency  of  the  A/D  clock  could  be  varied  to  optimize  the 
resolution  of  a  particular  frequency  range.  This  resolution  could  also  be  enhanced  by 
algorithms  which  include  skipping  and  multiple  FFTs.  Velocity  resolution  between  0.1% 
and  0.5%  should  thus  be  possible. 
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Recent  discussions  with  the  technical  program  monitor  indicate  that  measuring  the  particle 
size  would  considerably  enhance  the  capability  of  the  system.  It  is  our  opinion  and 
recommendation  that  such  capability  could  be  added  to  the  probe  without  a  substantial 
increase  in  its  complexity.  This  stems  from  the  fact  that  the  processor  is  programmable  and 
could  incorporate  additional  algorithms. 

2.  BACKGROUND  AND  THEORY 

The  LDV  transmitter  produced  three  laser  beams  which  crossed  to  form  the  probe  volume. 
Particles  traveling  though  this  probe  volume  scattered  light  which  contained  the  velocity 
information.  This  scattered  light  was  collected  by  a  receiver,  transformed  by  an  APD  into 
an  electronic  signal  with  similar  time  and  amplitude  variations,  and  digitized  and  stored  by 
the  electronic  processor.  To  measure  two  velocity  components  each  of  the  three  laser 
beams  was  frequency  shifted  by  a  Bragg  cell  such  that  the  difference  in  frequency  between 
any  two  beams  was  unique  and  separable  in  the  frequency  domain.  Figure  1  shows  the 
three  laser  beam  configuration  while  Figure  2  shows  a  schematic  of  the  breadboard  used  in 
Phase  I.  As  indicated  on  Figure  1,  the  three  frequencies  were:  75,  80,  and  90  MHz,  which 
produced  frequency  differences  of  5,  10,  and  15  MHz.  Thus  the  horizontal  velocity 
component  had  a  frequency  shift  (fs,)  of  5  MHz,  the  vertical  velocity  component  had  a 
frequency  shift  (fs2)  of  10  MHz,  and  tne  15  MHz  shift  corresponded  to  a  diagonal  velocity 
component  which  was  discarded.  A  stationary  particle  positioned  in  the  probe  volume 
would  scatter  light  which  contained  all  three  frequencies.  As  the  particles  moved  due  to 
flow  velocity,  the  processed  frequencies  would  be: 


fj-fcjiU  /Sp  (1) 

f2  =  fs2±V/52,  (2) 

where  f,  and  f2  are  the  processed  frequencies,  5 2  and 6,  are  the  fringe  spacings,  and  U  and  V 
are  the  velocity  components  in  the  horizontal  and  vertical  directions.  The  values  of  ^  and 
s ^  were  28  /im  and  25  /im  respectively.  U/5  is  referred  to  as  the  Doppler  frequency  and  its 
sign  can  be  easily  established  by  comparing  f  to  fs. 

The  above  theory  leads  to  the  conclusion  that  a  single  photodetector  and  a  single  electronic 
processor  can  provide  the  velocity  magnitude  and  direction  in  two  axes.  This  was 
successfully  demonstrated  with  the  Phase  I  system.  The  constraint  imposed  by  this  system 
is  that  f,  and  f2  must  be  separable  in  the  frequency  spectrum.  This  can  be  easily  met  in 
ocean  applications  since  the  velocities  are  very  slow.  A  very  conservative  design  would 
require  that  the  difference  between  fj  and  f2  be  four  times  the  maximum  expected  Doppler 
frequency. 

The  Phase  I  design  was  based  on  a  maximum  velocity  of  10  m/s  crossing  a  25  Mm  fringe 
spacing,  therefore  yielding  a  Doppler  frequency  of  0.4  MHz.  Thus,  the  difference  between 
processed  frequencies,  must  be  fj  -  f  =  1.6  MHz.  This  can  be  achieved  with  Bragg  cells 
operating  at  80  MHz,  81.6  MHz,  and  84.8  MHz.  The  maximum  processed  frequency  would 
then  be  5.2  MHz  (for  a  0.4  MHz  shift)  and  the  minimum  frequency  would  be  1.2  MHz.  A 
20  MHz  digitizer  would  exceed  the  requirements  imposed  by  Nyquist  limit,  and  provide 
excellent  velocity  resolution.  During  Phase  I  we  used  a  100  MHz  digitizer  which  is  much 
faster  than  necessary  and  uses  more  power  than  a  20  MHz  digitizer,  which  according  to  this 
simple  analysis  would  be  more  adequate  for  the  Phase  II  system. 

For  most  ocean  applications  the  velocity  would  be  even  slower  than  the  value  assumed 
above;  probably  in  the  order  of  cm/s.  In  this  case  the  digitizer  could  be  much  slower  than 
the  chosen  20  MHz  resulting  in  even  lower  power  consumption. 
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'inure  2.  Phase  I  breadboard. 


Mie  scattering  solutions 


The  scattering  cross-section  of  particles  between  approximately  2  and  20  microns  in 
diameter  were  computed  for  three  indices  of  refraction  and  two  collection  angles.  The 
indices  of  refraction  were  1.05,  1.1,  and  1.2  which  according  to  the  literature  are  fairly 
typical  of  particles  found  in  the  ocean.  The  angles  of  collection  were  160  and  180  degrees. 
While  the  Phase  II  system  would  operate  at  180  degrees,  the  Phase  I  collection  angle  was 
chosen  as  160.  This  was  due  to  the  fact  that  we  were  short  of  one  beam  splitter  to 
implement  the  180  degree  configuration.  For  the  purposes  of  system  evaluation  it  makes 
little  difference  since  the  scattering  cross-section  is  or  the  same  order  magnitude  for  the 
two  collection  angles.  The  other  parameters  which  influence  the  scattering  cross-section 
are  the  laser  wavelength,  the  F/#  of  the  receiving  lens,  and  the  polarization  of  the  incident 
beam.  These  parameters  are  listed  with  the  results. 

The  results  are  given  on  Tables  1  through  6  and  on  Figures  3  and  4.  These  results  were 
used  in  predicting  signal  to  noise  ratios.  It  is  evident  from  these  figures,  that  the  scattered 
light  intensity  collected  in  back  scattering  is  a  strong  function  of  index  of  refraction.  This 
combined  with  the  nonpherical  shape  of  the  particles  precludes  accurate  particle  sizing 
using  backscatting.  An  on-axis  configuration  would  collect  light  which  is  scattered  by 
diffraction  which  is  more  insensitive  to  shape  and  index  of  refraction  and  would  therefore 
yield  more  accurate  results.  Appendix  A  and  B  show  possible  schemes  which  would 
provide  the  particle  velocity  and  size  distributions.  It  is  shown  that  for  angles  of  collection 
near  on-axis  (zero  degrees)  the  intensity  of  the  scattered  light  is  reasonably  independent  of 
index  of  refraction. 
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DATA  FILENAME  :onl20160 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RADIUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  160.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.200  .000 


F (1)  :  RESPONSE  FUNCTION 

I/IO  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 
CROSS  SECTION  (1.0E-8  CMA2) 


SIZE  PARA. 
10.0000 
20.0000 
30.0000 
40.0000 
50.0000 
60.0000 
70.0000 
80.0000 


DIAM. (MICRON) 
2.5942 
5.1885 
7.7827 
10.3769 
12.9711 
15.5654 
18.1596 
20.7538 


F  ( 1) 

. 14 15E+00 
. 2661E+01 
. 5932E+01 
. 8399E+01 
. 2166E+02 
. 1147E+02 
. 3338E+02 
. 2816E+02 


I/IO 

. 2380E-02 
. 4477E-01 
. 9981E-01 
. 1413E+00 
. 3644E+00 
. 1930E+00 
. 5617E+00 
. 4739E+00 


Table  1.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.2 
collected  at  160  degrees. 


DATA  FILENAME  :onll0160 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RAPTUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  160.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.100  .000 


F (1)  :  RESPONSE  FUNCTION 

I/IO  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 
CROSS  SECTION  (1.0E-8  CMA2) 


SIZE  PARA. 
10.0000 
20.0000 
30.0000 
40.0000 
50.0000 
60.0000 
70.0000 
80.0000 


DIAM. (MICRON) 
2.5942 
5.1885 
7.7827 
10.3769 
12.9711 
15.5654 
18.1596 
20.7538 


F  ( 1) 

. 3949E-01 
. 9905E-01 
. 2  639E+00 
. 9807E+00 
. 2854E+01 
.672  5E+01 
. 1181E+02 
. 2  048E+02 


I/IO 

. 6644E-03 
. 1666E-02 
. 4441E-02 
. 1650E-01 
. 4802E-01 
. 1131E+00 
. 1988E+00 
. 3447E+00 


Table  2.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.1 
collected  at  160  degrees. 
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DATA  FILENAME  :onl05160 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RADIUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  160.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.050  .000 


F ( 1)  :  RESPONSE  FUNCTION 

I/I 0  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 
CROSS  SECTION  (1.0E-8  CMA2) 


SIZE  PARA. 
10.0000 
20.0000 
30.0000 
40.0000 
50.0000 
60.0000 
70.0000 
80.0000 


DIAM. (MICRON) 
2.5942 
5.1885 
7.7827 
10.3769 
12.9711 
15.5654 
18.1596 
20.7538 


F  ( 1) 

. 9050E-03 
. 2273E-01 
. 3 192E-01 
. 5256E-01 
. 1078E+00 
. 1442E+00 
. 1878E+00 
. 3363E+00 


I/IO 

. 1523E-04 
. 3824E-03 
. 5371E-03 
. 8843E-03 
. 1814E-02 
. 2426E-02 
. 3160E-02 
. 5658E-02 


Table  3.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.05 
collected  at  160  degrees. 


DATA  FILENAME  :onl20180 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RADIUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  180.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.200  .000 


F ( 1)  :  RESPONSE  FUNCTION 

I/IO  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 
CROSS  SECTION  (1.0E-8  CMA2) 


SIZE  PARA. 
10.0000 
20.0000 
30.0000 
40.0000 
50.0000 
60.0000 
70.0000 
80.0000 


DIAM. (MICRON) 
2.5942 
5.1885 
7.7827 
10.3769 
12.9711 
15.5654 
18.1596 
20.7538 


F  ( 1) 

. 1146E+00 
. 7788E+01 
. 1741E+02 
. 3226E+02 
. 484  3E+02 
. 3 197E+02 
. 1245E+03 
. 1178E+03 


I/IO 

. 1928E-02 
. 13 lOE+OO 
. 2929E+00 
. 5427E+00 
. 8 14  8E+00 
. 5380E+00 
. 2094E+01 
. 1982E+01 


Table  4.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.2 
collected  at  180  degrees. 
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DATA  FILENAME  :onll0180 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RADIUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  180.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.100  .000 


F ( 1)  :  RESPONSE  FUNCTION 

I/IO  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 
CROSS  SECTION  (1.0E-8  CMA2) 


SIZE  PARA. 
10.0000 
20.0000 
30.0000 
40.0000 
50.0000 
60.0000 
70.0000 
80.0000 


DIAM. (MICRON) 
2.5942 
5.1885 
7.7827 
10.3769 
12.9711 
15.5654 
18.1596 
20.7538 


F(l) 

. 54  57E-01 
. 2808E+00 
. 9971E+00 
. 3010E+01 
. 7700E+01 
. 1744E+02 
. 3293E+02 
. 5105E+02 


I/IO 

. 9182E-03 
. 4724E-02 
. 1678E-01 
. 5065E-01 
. 1296E+00 
. 2935E+00 
. 554 1E+00 
. 8589E+00 


Table  5.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.1 
collected  at  180  degrees. 


DATA  FILENAME  :onl05180 

LENS  DIAMETER  IN  (MM)  =  20.000 

NUMBER  OF  INTEGRATION  ALONG  RADIUS  &  PI  =  5  5 

LENS  CENTER  DISTANCE  IN  (MM)  =  50.000 

LENS  CENTER  ANGLE  IN  (DEGREE)  =  180.000 

LASER  WAVELENGTH  IN  (MICRON)  =  .8150 

LASER  POLA  TO  X  (X  TO  Y)  IN  (DEGREE) =  90.000 

INDEX  OF  REFRACTION  =  1.050  .000 


F ( 1 )  :  RESPONSE  FUNCTION 

I/IO  :  RATIO  OF  SCATTER  TO  INCIDENT  INTENSITY 


CROSS 

;  SECTION  (1.0E- 

•8  CMA  2 ) 

SIZE  PARA. 

DIAM. (MICRON) 

F  ( 1 ) 

I/IO 

10.0000 

2.5942 

. 2668E-02 

.  4490E-04 

20.0000 

5.1885 

. 7440E-02 

. 1252E-03 

30.0000 

7.7827 

.  7256E-01 

.  1221E-02 

40.0000 

10.3769 

. 392  IE-01 

. 6597E-03 

50.0000 

12.9711 

. 6340E-01 

. 1067E-02 

60.0000 

15.5654 

. 3149E+00 

. 5298E-02 

70.0000 

18.1596 

. 2 17 1E+00 

.  3652E-02 

80.0000 

20.7538 

. 3084E+00 

. 5189E-02 

Table  6.  Mie  scattering  computations  for  particles  with  an  index  of  refraction  of  1.05 
collected  at  180  degrees. 
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Scattering  vs  Size  Parameter  at  a  160°  Collection  Angle 


1 .000 


0.100 


0.010 


-  0.001 


1 .000E-4 


1.000E-5 

0  20  40  60  80  1C 

Size  Parameter  (dimensionless) 


Figure  3.  Mie  scattering  computations  for  particles  with  three  indices  of  refraction 
collected  at  160  degrees. 
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3.  EXPERIMENTAL  BREADBOARD 

The  key  elements  of  the  Phase  I  breadboard  (see  Figure  2)  are: 

*  Laser  diode  Melles  Griot  model  06DLL707  with  temperature  compensation  and  power 
control. 

*  Three  Bragg  cells  Isomet  model  1205C-1  specially  coated  to  cover  the  laser  diode 
wavelength  of  815  nm.  The  Bragg  cells  were  powered  by  three  separate  drivers  with 
corresponding  constant  frequencies  of  75,  80,  and  90  MHz.  This  choice  of  frequencies 
allowed  for  spectral  separation  of  both  the  horizontal  and  vertical  velocity  components. 

*  Avalanche  photodiode  (APD)  Texas  Optoelectronics  model  TIED88  and  custom  made 
preamplifier  board. 

*  Programmable  digital  signal  processor  (DSP). 

*  Optics  to  produce  probe  volume  and  collect  the  scattered  light. 

The  laser  diode  operated  TEM^at  815  nm  and  was  capable  of  providing  40  mw  of  power. 
For  the  purposes  of  Phase  I  the  laser  diode  included  a  commercial  temperature  controller 
and  power  regulator  which  are  bulky  and  not  appropriate  for  an  underwater  system. 
However,  these  controllers  would  not  be  necessary  in  an  ocean  environment  since  the 
temperature  of  the  ocean  does  not  change  by  more  than  about  20  F  in  any  particular 
location. 

The  laser  beam  was  split  into  three  beams  and  each  was  frequency  shifted  by  one  of  the 
Bragg  cells.  These  beams  were  crossed  by  a  transmitting  lens  to  form  three  overlapping 
fringe  patterns  which  were  spectrally  separated.  The  difference  in  frequency  between  the 
two  horizontal  beams  and  the  two  vertical  beams  was  5  MHz  and  10  MHz  respective^.  A 
parasitic  diagonal  shift  of  15  MHz  is  also  obtained  but  it  was  easily  filtered  out. 

The  APD  is  the  optimum  photodetector  for  near  infrared  detection.  Its  power  requirement 
is  much  smaller  than  that  required  by  a  photomultiplier  tube  since  its  quantum  efficiency  is 
about  ten  times  larger.  Furthermore  the  APD  is  small  and  rugged.  Unlike  conventional 
laser  Doppler  velocimeters  in  which  one  detector  is  needed  for  each  velocity  component,  a 
single  detector  was  employed  here  to  measure  the  vertical  and  horizontal  velocity 
components.  This  is  possible  because  the  signals  can  be  spectrally  separated  as  explained 
above. 

The  programmable  electronic  processor  was  made  by  Signatec  and  included  a  Texas 
Instrument  DSP320C25  which  capable  of  107  operations  per  second.  Thus  fast  Fourier 
transforms  can  be  performed  in  a  few  milliseconds.  Furthermore,  the  processor  could  be 
programmed  to  perform  any  combination  of  algorithms  for  signal  acquisition  and 
processing.  While  the  Phase  I  processor  was  hosted  by  an  IBM  personal  computer,  the 
Phase  II  version  would  stand  alone. 

The  laser  beam,  which  was  nominally  8  mm  in  diameter,  was  minified  to  about  1  mm  by  a 
combination  of  two  lenses  (L,  and  L,)  with  focal  lengths  of  200  and  25.4  mm.  The  beam 
splitter  BS,  separated  the  beam  into  two  parts,  a  reflected  part  and  a  transmitted  part.  The 
transmitted  portion  of  the  beam  was  split  again  by  BSr  Each  one  of  these  three  beams 
entered  the  respective  Bragg  cell.  A  combination  of  mirrors  and  lens  L3  was  used  to  cross 
and  focus  the  beams  into  a  probe  volume  of  about  60  nm  in  diameter  and  with  vertical 
fringe  spacing  of  about  25  /im,  and  horizontal  fringe  spacing  of  28  ^m,  The  light  scattered 
by  particles  was  collected  at  160  degrees  and  focused  on  the  APD  by  two  lenses  L4  and  L5) 
with  an  equivalent  F/2.5.  The  APD  converted  the  scattered  light  into  an  electronic  signal 
which  was  processed. 
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To  demonstrate  the  three  beam  laser  Doppler  velocimeter  in  water,  a  Plexiglass  test 
chamber  was  constructed.  A  pump  flowed  water  through  the  chamber  and  a  2  mm  orifice 
and  valve  were  used  to  control  the  flow  velocity.  With  the  valve  wide  open  the  flow  rate 
was  measured  to  be  1.5  liters/min,  corresponding  to  an  average  velocity  at  the  2  mm  orifice 
of  approximately  8  m/sec.  Polystyrene  particles  with  a  diameter  of  9.87  /im  and 
monodisperse  to  within  0.057 nm  were  added  to  the  water  to  provide  calibrated  scatterers 
to  be  measured  by  the  LDV. 

4.  RESULTS 

Signal  amplitude  and  signal  to  noise  ratio 

The  largest  signals  corresponded  to  particles  traveling  through  the  center  of  the  Gaussian 
probe  volume.  The  peak  pedestal  of  these  signals  had  values  of  about  200  mv.  The 
amplitude  of  these  signals  could  be  predicted  with  the  following  expression: 

Signal  Amplitude  =  lpxR, 

where  I  represents  the  peak  of  the  pedestal  and  is  given  in  watts  of  collected  light,  and  R  is 
the  detector  response  which  is  given  in  volts/watt. 

I  was  estimated  from  the  Mie  scattering  solution  knowing  the  particle  size  (10  /im),  the 
ifidex  of  refraction  (1.2),  the  angle  of  collection  (160),  the  F#  or  collection  (2.5),  and  the 
laser  wavelength  (815  nm).  The  scattering  cross-section  of  a  10  micron  polystyrene  particle 
in  water  was  thus  computed  to  be  about  0.14  /im2  (see  the  value  of  I/Io  on  Table  1  for  a 
10.37  /*m  particle).  TTie  illuminating  beams  had  a  combined  total  power  of  20  mw  and 
focused  to  a  60  /im  diameter  probe  volume.  The  power  density  in  the  middle  of  the  probe 
volume  was  then  estimated  at  7  x  10-6  w/(im2.  The  power  density  multiplied  by  the 
scattering  cross-section  yields  the  total  collected  light,  which  is  9.8  x  10  '  watts. 

The  response  of  the  detector  is  given  by: 

R  =  Re  x  APD  gain  x  terminal  resistance  x  amplifier  gain. 

For  the  APD  at  a  wavelength  of  815  nm,  the  radiant  responsivity  R  =  0.3  amps/watt;  the 
APD  gain  was  calibrated  to  be  25,  the  terminal  resistance  was  1000  onms,  and  the  amplifier 
gain  was  75.  These  values  yield  a  detector  response  of  5.6  x  105  volts/watt. 

The  peak  signal  amplitude  can  then  be  computed  to  be  550  mv.  The  signal  output  was  ac 
coupled  which  limited  the  peak  signal  to  about  50%  of  its  value,  that  is  275  mv.  This 
prediction  is  very  close  to  the  measured  value  of  200  mv  and  gave  us  confidence  to  predict 
the  signal  level  associated  with  particles  of  different  diameter  and  index  of  refraction. 

For  the  amplifier  gain  of  75  the  peak  to  peak  electronic  noise  was  about  40  mv.  This  was 
established  given  that  an  ac  threshold  of  20  mv  excluded  triggering  on  the  electronic  noise. 
Since  the  measured  signals  had  a  peak  amplitude  of  200  the  signal  to  noise  ratio  was 
about  10.  This  signal  to  noise  ratio  could  also  be  established  from  the  frequency  spectrum 
by  comparing  the  height  of  the  signal  bins  and  the  adjacent  noise  bins.  It  was  desirable  to 
set  the  threshold  level  above  the  electronic  noise  to  avoid  triggering  on  the  noise  and  fill  up 
the  memory  with  useless  information.  Since  this  signal/noise  corresponded  to  a  scattering 
cross  section  of  0.14  /an2,  it  could  be  concluded  that  the  minimum  scattering  cross-section 
would  be  approximately  0.014  /im2  for  a  signal  to  noise  ratio  of  unity  This  minimum  could 
be  decreased  by  reducing  the  noise,  and  there  are  various  ways  to  meet  this  goal.  First, 
since  the  frequency  range  of  interest  was  established  to  be  between  1.1  and  5.5  MHz, 
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analog  filters  could  be  used  to  limit  the  frequency  to  this  range.  This  would  reduce  the 
bandwidth  to  about  10  MHz  instead  of  the  50  MHz  used  during  Phase  I.  Second,  more 
efficient  APD  and  amplifier  combinations  could  be  selected.  The  combined  effect  of  these 
two  improvements  would  result  in  reducing  the  noise  by  about  a  factor  of  10.  Under  these 
conditions  the  minimum  scattering  cross-section  would  be  0.0014  /im2. 

The  scattering  cross-section  is  a  function  of  particle  size  and  index  of  refraction.  Tables  1 
through  6  show  sample  computations  obtained  by  solving  the  Mie  scattering  equations  for 
indices  of  refraction  1.05,  1.1,  and  1.2.  These  values  are  typical  of  the  micro-organisms  and 
sediment  found  in  the  ocean.  If  we  use  the  value  of  0.0014  nm2  as  the  minimum  scattering 
cross  section  detectable  by  the  system,  it  can  be  seen  from  these  tables  that  for  an  index  of 
refraction  of  1.2  the  minimum  detectable  diameter  is  about  2  ^m;  for  an  index  of  refraction 
of  1.1  the  minimum  detectable  diameter  is  about  5  jim,  and  for  an  index  of  refraction  of 
1.05  the  minimum  detectable  diameter  is  about  12 /im. 

The  velocity  distribution 

Experiments  were  run  at  various  flow  rates  and  with  the  test  chamber  inclined  at  various 
angles  to  produce  different  velocity  components  in  the  horizontal  and  vertical  directions. 
Figures  5  and  6  show  the  results  of  a  sample  burst  corresponding  to  one  of  these 
experiments  for  an  inclination  of  the  test  chamber  of  34  degrees  and  the  flow  control  valve 
wide  open.  The  frequency  spectrum  was  obtained  by  performing  an  FFT  on  the  collected 
data.  Notice  the  values  of  fj  =  5.225  MHz  and  f2  =  10.25  MHz.  Substituting  these  values 
into  equations  (1)  and  (2),  we  arrive  with  velocity  values  of  U  =  V  =  6.2  m/s.  In  this 
particular  example  the  U  and  V  components  are  equal  and  correspond  to  flow  moving  at 
45  degrees.  In  general,  U  and  V  can  have  arbitrary  values.  This  arrangement  permits  the 
measurement  of  flow  velocity  and  turbulence  at  arbitrary  directions;  including  recirculating 
flows.  Figure  7  shows  a  portion  of  the  signal  corresponding  to  a  stationary  particle  (pump 
off).  The  Fourier  transform  of  this  signal  is  shown  on  Figure  8  and  demonstrates  that  for 
the  no  flow  case  the  measured  frequency  corresponds  only  to  the  Bragg  cell  frequency  shift 
(i.e.  5  and  10  MHz  respectively). 

A  sample  of  13  bursts  was  analyze  to  study  the  turbulence  of  the  flow  coming  out  of  the 
orifice.  This  is  a  very  small  sample  for  any  statistical  purposes  but  it  serves  to  demonstrate 
a  concept.  The  mean  velocity  of  the  flow  at  this  point  based  on  the  13  measurements  was 
8.76  m/sec.  The  direction  of  flow  was  found  to  be  56°.  The  standard  deviation  of  both  the 
velocity  and  direction  of  the  flow  was  found  to  be  10%,  indicating  the  high  degree  of 
turbulence.  A  visual  display  of  these  measurements  is  shown  on  Figure  9. 

PHASE  II  SYSTEM 

For  a  maximum  velocity  of  10  m/s  and  assuming  a  fringe  spacing  of  25  microns  we 
estimated  above  that  the  required  frequency  shifts  would  be  1.6,  and  3.2  MHz.  These  could 
be  achieved  with  three  Bragg  cells  operating  at  80,  81.6,  and  84.8  MHz,  or  with  two  Bragg 
cells  operating  at  80  and  81.6  MHz  while  utilizing  the  first  and  second  orders. 
Alternatively,  these  frequency  shifts  could  be  obtained  by  running  three  laser  diodes  at 
slightly  different  wavelengths.  In  Phase  I  we  demonstrated  that  the  three  Bragg  cell 
approach  is  sound;  the  two  Bragg  cell  approach  appears  to  be  a  straight  forward  extension. 
Here,  two  Bragg  cells  would  be  driven  with  frequencies  of  80  MHz  and  81.6  MHz,  thus  the 
second  order  frequencies  would  be  160  and  163.2  MHz.  The  horizontal  velocity 
component  could  be  measured  by  mixing  the  first  order  beams  which  yield  a  frequency  shift 
of  1.6  MHz;  the  vertical  component  of  the  velocity  could  be  measured  by  mixing  the  second 
order  beams  which  yield  a  frequency  shift  of  3.2  MHz.  It  may  also  be  worth  considering  the 
three  laser  approach  in  the  beginning  of  Phase  II  since  if  they  prove  successful,  they 
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consume  less  power  than  the  Bragg  cells.  Since  this  approach  has  not  been  demonstrated, 
we  will  describe  the  Phase  II  system  under  the  assumption  that  three  Bragg  cells  will  be 
used  for  frequency  shifting. 

The  complete  Phase  II  system  is  expected  to  measure  7"  x  26"  as  indicated  on  Figure  10. 
The  system  will  include  a  battery  pack,  the  electronic  processor,  memory,  laser  power 
supply  (this  may  or  may  not  be  necessary),  Bragg  cell  drivers,  single  mode  diode  laser, 
telescope  for  beam  conditioning,  two  beam  splitters,  three  Bragg  cells,  three  rhombs  for 
beam  steering,  a  plate  beam  splitter,  transmitting  lens,  receiving  lens,  and  avalanche 
photodiode.  The  shape  of  the  system  is  expected  to  be  cylindrical  with  a  single  window  for 
transmitting  the  laser  beams  and  receiving  the  backscattered  light. 

An  alternative  Phase  II  system  is  shown  on  Figure  11.  This  alternative  system  would 
incorporate  particle  sizing  capability.  Based  on  the  discussion  above,  the  beast  result  are 
expected  when  the  light  is  collected  in  the  forward  scattering  direction.  Thus,  this  system 
would  include  a  fiber  optic  receiver  which  collects  the  scattered  light  and  transfers  it  to  the 
photodetector.  Furthermore,  the  probe  volume  may  have  to  be  shaped  as  an  elongated 
sheet  as  indicated  on  Appendix  A,  to  avoid  the  Gaussian  ambiguity  associated  with  a  round 
beam.  The  final  configuration  will  depend  on  the  needed  resolution.  The  alternative 
system  will  also  include  the  capability  of  measuring  the  transmitted  laser  beam  in  order  to 
correct  the  scattered  signal  by  the  obscuration. 

Data  acquisition  and  memory  requirements. 

The  choice  between  solid  state  memory  and  a  more  permanent  record,  such  as  a  tape, 
would  have  to  be  made  during  the  Phase  II  effort.  Solid  state  memory  requires  very  little 
power  to  operate  but  it  is  vulnerable  to  power  outage.  A  tape  backup,  or  some  other 
means  of  permanent  memory  may  be  necessary  if  the  chances  of  power  interruption  are 
reasonably  large.  The  information  to  be  stored  in  memory  consists  of  the  mean  and 
fluctuation  of  the  horizontal  and  vertical  velocity  components,  as  well  as  the  Reynolds 
stress  and  the  time  of  collection.  A  total  of  six  2-byte  parameters  (U  bar,  V  bar,  U\  V’, 
U’V\  t)  are  expected  to  be  stored  for  each  run.  A  run  should  consist  of  approximately 
10,000  particles  measured  over  20  seconds  or  so  every  hour  or  as  defined  by  a  data 
collection  logic.  The  time  of  collection  or  the  data  rate  would  have  to  be  automatically 
adjusted  to  account  for  variations  on  the  data  rate.  The  raw  data  associated  with  10,000 
particles  is  about  1  Mb.  To  avoid  a  large  electronic  processor  actual  high  speed  memory 
will  be  limited  to  64  Kb  and  the  data  will  be  acquired  in  16  intervals.  The  information 
associated  with  these  10,000  particles  will  then  be  condensed  to  six  2-byte  parameters  for  a 
total  of  12  bytes. 

Assuming  that  data  are  acquired  every  hour  for  six  months,  a  total  of  4320  runs  will  be 
stored  which  translates  to  51840  bytes.  Alternatively  data  may  be  acquired  when  prompted 
by  an  external  sensor  such  as  a  microphone  or  an  accelerometer.  Under  this  last  scenario 
data  would  be  acquired  when  the  conditions  of  the  ocean  exceed  some  pre-determined 
threshold. 
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Drivers 


Alternative  Phase  II  LDV  system  utilizing  forward  scattering 


6  watts 


Power  budget 

20  MHz  electronic  processor  with  64  Kb  of  memory 


Three  Bragg  cells  3  watts 

40  mw  laser  diode  0.2  watts 

APD  0.16  watts 

External  sensor  (if  chosen)  0.02  watts 


The  total  power  consumption  would  be  about  10  watts  during  data  acquisition,  and  about 
0.02  watts  during  waiting  periods. 

For  24  events  per  day  with  a  duration  of  20  seconds  the  total  data  acquisition  time  for  six 
months  would  be  24  hours.  The  waiting  time  would  be  4296  hours.  The  total  energy 
required  by  the  system  would  then  be: 

24  hrs  x  10  watts  +  4296  hours  x  0.02  watts  =  326  watt-hours. 

The  system  could  be  designed  to  operate  at  12  volts  thus  requiring  27  amp-hours.  There 
are  several  batteries  which  appear  to  satisfy  this  energy  requirement  while  maintaining  a 
reasonable  size.  Examples  of  these  batteries  include  lead  calcium  and  sealed  lead  acid. 
The  precise  choice  of  batteries  will  be  an  important  task  during  Phase  II. 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  The  programmable  processor  allows  the  implementation  of  intelligent  algorithms 
which  can  be  specifically  tailored  to  optimize  the  data  acquisition  and  processing  in  a 
marine  environment. 

2.  The  three  Bragg  cell  concept  allows  for  simple  discrimination  of  two  velocity 
components  using  a  single  wavelength  laser  and  a  single  photodetector. 

3.  The  diode  laser  provides  ample  power  with  very  little  energy  consumption  making  it  an 
ideal  choice  for  a  battery  operated  system. 

4.  Particle  sizing  can  be  added  to  the  Phase  I  system  with  reasonably  simple 
modifications.  Forward  scattering  collection  may  be  required  to  obtain  adequate 
signal  resolution. 

To  demonstrate  the  three  beam  LDV  in  water,  a  Plexiglass  test  chamber  was  constructed. 
A  pump  flowed  water  through  the  chamber  and  a  2  mm  orifice  and  valve  were  used  to 
control  the  flow  velocity.  With  the  valve  wide  open  the  flow  rate  was  measured  to  be  1.5 
liters/min,  corresponding  to  an  average  velocity  at  the  2  mm  orifice  of  approximately  8 
m/sec.  Polystyrene  particles  with  a  diameter  of  9.87 /tm  and  monodisperse  to  within  0.057 
Atm  were  added  to  the  water  to  provide  calibrated  scatterers  to  be  measured  by  the  LDV. 
Two  lenses  were  used  to  provide  an  F/2.5  collection  system  which  focused  the  scattered 
light  onto  an  APD  detector.  The  signal  from  this  detector  was  amplified  and  sent  to  the 
programmable  DSP. 
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Abstract 

A  technique  combining  the  visibility  of  a 
Doppler  signal  and  the  intensity  of  the  scattered 
light  to  measure  the  size  and  velocity  of  parti¬ 
cles  is  presented.  It  is  shown  that  using  only 
the  visibility  technique  can  lead  to  large  errors 
under  many  conditions  such  as  dense  sprays.  It  is 
also  shown  that  this  error  is  considerably  reduced 
and  very  high  resolution  is  obtained  by  combining 
the  visibility  with  the  intensity  of  the  scattered 
light.  An  instrument  was  developed  using  this  new 
concept  and  measurements  were  performed  in  sprays 
of  known  .Laracteri  sLics,  The  results  of  monod is- 
persed,  hi  modal  and  trimodal  sprays  are  reported. 

I  nt  roduct ion 

The  need  for  characterizing  droplet  environ¬ 
ments  (especially  sprays)  in  many  combustion  pro¬ 
cesses  is  well  established.  Accurate  spatially 
resolved  measurements  of  the  size,  velocity  and 
number  density  of  the  droplets,  and  the  velocity 
of  the  gas  are  needed  by  system  developers  and 
mode  lers . 

The  ability  to  distinguish  between  droplet 
velocity  and  gas  velocity  is  important  since 
together  with  the  droplet  size  it  yields  local 
Reynolds  numbers  from  which  particle  drag  coeffi¬ 
cients  can  be  determined. 

In  a  combustion  environment,  the  rate  of  evap¬ 
oration  of  the  droplets  is  determined  by  the  rate 
of  change  of  diameter  of  individual  droplets. 
This  quantitative  information  of  the  drop  size  and 
velocity  distributions  is  needed  throughout  the 
spray  as  a  function  of  space  and  time. 

Techniques  are,  therefore,  needed  that  can 
provide  the  above  information  noni nt rus i ve 1 v , 
accurately  and  in  real  time  from  real  sprays. 
There  are  several  techniques  that  meet  some  of  tin- 
conditions  listed  above.  They  can  be  grouped  as 
imaging  and  scattering  techniques.  However,  only 
a  few  of  these  techniques  can  be  used  for  simul¬ 
taneous  measurement  of  size  and  velocity  at  a 
point.  The  two  parameters  used  to  date  in  these 
techniques  are  the  visibility  ol  a  Doppler  signal 
and  the  absolute  intensity  of  the  scattered  light. 
This  paper  reports  very  recent  research  that  shows 
the  errors  associated  with  the  visibility  tech¬ 
nique,  and  a  new  method  that  uses  simultaneously 
the  visibility  ot  the  Doppler  signal  and  the  pede¬ 
stal’s  peak  intensity  to  perform  a  size  measure¬ 
ment,  This  technique  is  referred  t  o  is 

Visibility/Intensity  (V/I). 

Results  obtained  with  the  V/I  technique  . in- 
reported  for  monod  i  spe  rsed ,  himndul  and  trir.mdal 
sprays  and  monod i spe rsed  string  ot  droplets.  The 
mond i spe rsed  sprays  offer  great  potential  as  cali¬ 
bration  sources.  The  results  show  that  hv  virtue 
of  <•  ombi  ui  :i>*  the  visibility  and  intensity  v»-rv 
precise  measurements  are  possible  in  sprav*.  »ud 

toforijjhl  \mrri*  ;in  IrvUihilr  ol  \<-r*Hi,iufi< ■«.  and 
\*Urnn;»iiht lm  .  I‘tt4.  \H  ni»ho  ri  M-n ul 


environments  where  visibility  only  -v  i  1 1  not 
perform. 

Associated  Research 

The  i  nter  f  erorW  ri  c  technique  to  measure  tin* 
size  and  velocity  ot  particle.*  flowing  in  a  fluid 
has  been  in  existence  tor  approximately  one 
decade.  Farmer  introduced  the  concept  of  using, 
visibility  for  particle  sizing  in  the  forward 
scattering  direction.  Since  then,  many 

researchers  have  contributed  to  the  development  of 
the  technique.  Robinson  and  Clm  ,  for  instance, 
rederived  Farmer’s  results  using  a  more  rigorous 
approach,  and  the  definition,  of  visibility  through 
the  first  order  Bessel  function  w > *  confirmed 
experimentally  in  forward  scatter.  However,  con¬ 
trary  to  what  was  thought  at  the  time,  this  visi¬ 
bility  relationship  could  not  be  used  at  any  angle- 
other  than  the  forward  or  zero  degree.  Adrian  and 
Or  lot f 5  showed  that  the  simple  relationship  valid 
in  forward  scattering  could  not  be  used  in  back- 
scattering.  This  finding  was  also  confirmed  by 
Roberds  ,  who  in  addition  showed  that  the  config¬ 
uration  of  the  receiving  optics  would  also  affect 
the  relationship  between  visibility  and  size. 

The  biggest  known  limitation  of  the  visibility 
technique  was  its  inability  to  measure  particles 
in  a  dense  field.  The  main  reason  being  that  the 
probe  volume  when  observed  in  the  forward  direc¬ 
tion  is  big  and,  therefore,  the  criterion  of  sin¬ 
gle  particle  measurment  is  not  always  mot.  In 
addition,  the  technique  dictated  that  the  largest 
measurable  droplet  was  of  the  order  of  the  fringe 
spacing.  This  imposes  a  limitation  of  about 
200  um  to  the  largest  measurable  particle  and  also 
results  in  a  large  probe  volume,  since  a  mini  nun 
number  of  fringes  is  needed  (typically  8)  to 
accurately  process  a  Doppler  signal. 

Bat-halo5  formulated  a  mathematical  model  based 
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the  L  *mu-  es t  ab l  is hes  .i  i,n  ot  the  big  oiii* 

whem  tlu*  hyperbolic  viriition  of  Lh«-  visibility 
is  iiv  g  L  i  g  i  b  l  e .  Also  digital  mi  lysis  of  tlu.*  sig- 
•mIs  he  l  ps  establish  son**  t  luit  are  i  a  error.  The 
rose.  1 1  v  ii  report  ed  iier«*  shows  lli.it  ns  i  :iy  tiu*  inten¬ 
sity  of  the  jn'il*’st.iL  *>t  tlu*  scattered  light  in 
uMition  to  tin.*  vi-ibilitv  will  .*  1  i  -ii  nit  «*  many  of 
these  prints,  u  i  )»n  i  t  i  i*a  at  l  v  i  ::ip  r  ■  c  i  ng  measurement 
.  j  r .  ■  1 1  r  ■  1 1  •  v  . 

Krrors  in  tin-  Visibility  Measurement s 

Tlu-  pr»‘.i  i  v t  ion  or  tin*  droplet  size  based  on 
tin-  visibility  of  tlu*  light  scattered  by  the  drop- 
lot  iT'i.ssi:),:  an  i nt e it e i ome t r i c  pattern  of  fringes 
depends  on  tin*  pattern  itself.  Fer  a  perfect,  sys¬ 
tem  the  fringe  visibility  should  be  I'H)  and  the 
light  scattered  should  be  primarily  of  one  type 
(ret  rant  ion,  t  «*  f  Lee  t  i  on  or  d  i  1 t  tact  ion ) .  There 
are,  however,  many  causes  in  practical  environ¬ 
ments  ttuit  contribute  to  alter  the  fringe  pattern 
and,  therefore,  confuse  the  relationship  between 
visibility  and  size.  The  result  is  that  apparent 
si  :es  instead  of  true  sizes  are  often  measured. 
As  an  'sample,  droplets  interacting  with  the  laser 
bo  iris  bet  c.  re  they  cross  will  randomly  reduce  (or 
*  ,'en  destroy)  tae  fringe  visibility.  If  a  droplet 
is  "m  ms'  i  r**d  at  such  time,  it  will  appear  bigger 
than  its  true  size.  There  are  several  factors 
that.  j«rod'H*e  an  error  in  the  visibility  measure - 
•vnt  end  they  are  listed  below.  They  are  divided 
i  :t  5  two  *roups:  the  ones  that  reduce  tlu*  visi¬ 

bility  and  t. he  ones  that  increase  it. 

A)  !'  o' Lors  that  Reduce  the  Visibility 

1)  ^articles  prior  to  crossover  destroy 
fringe  contrast.  This  is  a  function 
of  the  spray  density. 

2)  Bean  excursions  due  to  turbulent  hot 
med  i  a . 

1)  Hyperholoida 1  reduction  of 
visibility. 

4 )  hi  gh  t  sea  1 1  e  red  by  re  f  lec  t  i  on / 

refraction  can  add  destructively  for 
la  rye  droplets  immersed  in  a  non- 
uniform  (typically  Gaussian)  beam. 

i)  'hilt  i  pie  particles  in  probe  volume. 

b )  Increase  in  signal  background. 

Normally  from  many  little  drops 
present  in  the  probe  volume. 

7)  Loss  of  fringe  contrast  due  to 

fiulty  components  (beamsplitter, 
etc.  ) 

B)  K.ictors  that  Increase  the  Visibility 

1)  Out  of  focus  drops  "masked”  by  the 
p i nhn le . 

.' )  High  time  rise  on  leading  edge  of 
s  igna  1  . 

Tlu*  Vi s  i  b i  1  i  t  y  /  I n t ens  i  ty  (V/I)  Technique 

This  technique  makes  use  not  only  of  the  visi¬ 
bility  of  a  Doppler  signal  but  also  the  peak 
intensity  of  the  pedestal.  Both  parameters  are 
available  in  the  signal  and  their  cross- 
correlation  can  be  used  to  eliminate  faulty 
signals  produced  in  many  practical  environments. 
This  technique  will  especially  prevent  small 


particles  (high  visibility)  l  roiu  up  pea  r  i  ng  as 
large.  Figure  l  shows  a  Doppler  t  r-io-  when  both 
the  visibility  and  the  peak  intensity  of  the 
pedestal  are  indicated.  Tlu*  re  is  a  relationship 
between  the  size  of  the  droplet  and  the  amount  of 
scattered  light  given  by  Mie  theory.  This  rela¬ 
tionship  can  he  used  to  eliminate  signals  with  an 
apparently  different  si/e  by  using  the  following 
logic:  droplets  that  produce  a  certain  visibility 

are  associated  with  a  given  size;  huice,  the” 
should  scatter  light  with  a  given  intensity 
(characterized  by  I  ).  Two  exceptions  are  contem¬ 
plated:  first,  droplets  with  the  correct  visibil¬ 

ity  will  scatter  different  amounts  of  Light  due  to 
t  he  Gaussian  nature  of  tlu*  probe  volume’s  inten¬ 
sity;  second,  droplets  with  .in  erroneous  visibil¬ 
ity  will  not  scatter  light  with  an  intensity 
corresponding  to  their  apparent  size. 

The  V/I  method  will  then  establish  intensity 
limits  for  (.-very  measured  visibility.  This  will 
produce  a  well  established  probe  volume  as  a  func¬ 
tion  of  size  and  will  reject  droplets  measured 
with  an  erroneous  visibility.  for  instance, 
assume  that  the  size  corresponding  to  a  visibility 

V,  is  d..  Then  the  light  scattered  will  have  a 
i  i  f  r  , 

pedestal  with  maximum  intensity  pj(max).  Since 

this  intensity  would  only  be  detected  from  drop¬ 
lets  crossing  through  the  middle  of  the  probe 
volume.  Figure  l,  more  relaxing  limits  are  pro¬ 
posed.  That  is  * p |  will  be  accepted  when 

^Pj(min)  •:  ^p  ^  p^(max).  obviously,  the 

broader  the  limits,  the  larger  the  error  that  can 
be  allowed  with  the  visibility  technique. 

Figure  2  shows  three  cases  that  illustrate  the 
correcting  properties  of  the  V/l  method.  Figure 
2a  shows  a  droplet  of  size  with  visibility  Vj 
and  intensity  p^.  Since  this  last  one  is  within 
established  limits,  the  signal  is  accepted  as 
valid.  Figure  2b  shows  the  signal  produced  by  a 
droplet  of  diameter  d,  but  the  visibility  was 
reduced  to  Vj  due  to  fringe  contrast  reduction. 
Since  the  visibility  is  V| ,  this  signal  would  be 
interpreted  as  having  a  diameter  d^.  However,  the 
peak  intensity  of  the  pedestal  p  corresponds  to 
its  true  size  (d,)  which  is  smaller  than  pj(min) 
(to  which  it  is  compared).  Therefore,  the  signal 
is  rejected.  Figure  2c  shows  the  signal  of 
another  droplet  of  diameter  d?  but  with  correct 
visibility  V,  and  corresponding  pedestal  p?. 
This  last  one  will  also  be  accepted. 

Since  the  intensity  limits  are  functions  of 
droplet  size,  a  variable  filter  like  the  one  shown 
on  Figure  3  is  necessary.  It  contains  the  inten¬ 
sity  limits  l  (min),  I  (max)  corresponding  to  any 
visibility  V.  Figures  4a  and  4b  show  the  raw  data 
and  software  filtered  data  of  an  i ntensi ty /visi¬ 
bility  plot.  Figure  *a  shows  that  the  intensity 
is  bounded  by  l  (max)  (with  a  few  exceptions  due 
to  light  reradiated  bv  upstream  droplets).  Points 
with  different  intensi ty /vi s i bi 1 i t y  combinations 
can  be  either  due  to  drop  lots  moving  through  the 
»*dge  of  the  Gaussian  intensity  profile  of  the 
probe  volume  or  to  droplets  with  an  erroneous 
visibility.  By  limiting  the  intensity  of  the 
pedestal  for  which  a  signal  is  considered  valid, 
two  things  are  accomplished:  (1)  tlu*  limits  of 

the  probe  volume  .ire  carefully  established;  and, 
(2)  most.  importantly,  most  of  the  signals  with 
erroneous  visibility  are  rejected. 


ReSU  Its 


Experimontai  result <  -v'e  presented  to  illus¬ 
trate  the  accuracy  and  resolution  or  the  visibil¬ 
ity 'intens  i  ty  technique,  and  ’now  it  compares  with 
visibility  only, 

Effect  nf  Bean  Blockade  on  Size  Dis*  rihut ion 


vibration  ut  the  orifice,  and  with  the  d  i  r  i  ■ 
air  a  spray  a  these  droplets  was  lamed.  Typical 

number  dens; tv  at  the  probe  volume  was  a1;-’  drop- 

1  ' 

lets -'em  ami  the  diameter  of  t:ie  spray  at  the 
:  lane  of  measurements  was  about  9  mm. 

A  US  I  system  with  VI  cap  ob  i  l  i  ty  was  used  to 
ibtaiu  tiie  data  shown  be  1  rw. 


The  effect  of  a  very  dense  spray  block!  m:  the 
laser  beams  before  the  probe  volume  was  studied. 
Figure  5a  shows  the  hi s  Centrum  produced  when  mono- 
dispersed  droplets  travel  thmujh  tin.*  middle  of  an 
undisturbed  probe  volume.  Figure  5b  shows  a  his¬ 
togram  obtained  with  the  visibility  technique  of 
the  same  nonod i spersed  droplets  but  now  the  lamr 
beams  have  traveled  throe, :L  a  real  spray  before 
•::r- >ss  ine .  Wo  made  sure  that  t he  spray  was  not 
■to  in-:  int)  the  probe  vulur/*.  The  frin.te  pattern 
pr  iduced  by  crossing  two  laser  beams  i  well 
behaved  and  t heoret  i  rally  predicted  in  the  uud’s- 
turbed  -ase,  However,  when,  the  beams  tra.ul 
thr  <u,;h  a  spray,  this  fringe  pattern  can  be  con¬ 
siderably  1 1 1*- red  due  to  random  beam  blockade. 

The  result  can  be  a  distorted  Doppler  siynai  with 
a  visibility,  both  lower  and  1 1 i  her  than  that 
predicted  theoretically.  Considerable  error  in 
the  size  distribution  will  then  bo  observed.  The 
size  of  the  err>r  depends  on  the  visibility 
itself.  It  is  quite*  apparent  from  the  histograms 
that  what  should  be  »  monod i spersed  distribution 
is  measured  as  a  very  broad  o.je.  The  next  histo¬ 
gram  shown  on  Fi cure  5.  shows  the  results  obtained 
with.  V  I  when  the  beams  travel  thr  njph  the  spray 
before  crossin  :  •  -same  conditions  as  Figure  5b). 

Notice  that  a  very  .narrow  distribution  is  thus 
obtained . 

Measurements  with  the  her.?  lurid  hi  u  or  op  let 

Gene  rat->r 

The  calibration  and  verification  source  us*  i 
in  these  experiments  was  the  ik  rplund  hiu  nonodis- 
perse  droplet  generator.  This  //.uier-itoi  <an  work 
in  two  modes:  fa)  it  can  generate  a  string  of 

droplets  of  known  size  and  equally  spaced  at  tie* 
measurement  point;  (',)  us i rip,  a  cap  and  dispersion 
air,  these  droplets  can  be  dispersed  into  a  spray 
'f  monod  i  s;  ,*rsed  dr. plots. 

Tt  has  been  observed  in  previous  studies*"  that 
under  some  dispersion  air  conditions,  some  >r  the 
d  r  j|»  1  *•  t  s  will  collide  and  torn  other  droplets  with 
double  or  even  triple  the  volume.  Tlierefore,  it 
the  nr  i  mar/  droplet  d  i  uja*  ■  r. .  ■  r  d  droplets  with 

diameter  d  (  and  5  !  *  an  be  produced.  It 

must  be  p'*i’:feri  out  tint  in  these  studies,  w<*  d  i  d 
n  n  t  obtain  ph- -I.  o/ r  a  ;>h  i  me.jsu  rer;j*not  s  to  verify  tin* 

presence  >5  Me-  S  *  ■  ■  '  >  t!  <1 .  l  V  V  droplets.  'lie  si  /•'  ol 

the  primary  implet.s  was  accurately  predicted  from 
the  flow  rite  and  freque-ov  of  vibration  of  the 
p i nho l e . 

On**  limit  if  ion  imposed  by  the  fter/lund  hiu  i  . 
that  it  cannot  produce  droplets  ;f  a  d7»*i,  s  i  /<  at 
arbitrary  spacing.  Si  tee  the  measurement  instru¬ 
ment  requires  the  prese-ue  of  only  one  draper  in 
the  probe  volume,  that  limits  the  size  droplet 
that  can  he  measured. 

The  p  r  i r  e  *  1 1 1  |  e  used  in  these  e/pefi:.  -fits  W-V.  to 
produce  a  strinp  of  lur,y*  monod i spersed  droplets 
to  f  a  1  i  hr  a  t  **  the  i  us  t  r  n«*  n  f  .  I  hen  smaller  d  I  ■  »p- 
1  ''!  '  >  Were  produced  hv  i  n<  pus  i  l\,\  tile  frequency  ,/f 


opt  ical  conf  i  v.'urat  i  o  ’. : 

•  Co 1 1 v  c  t i on  a n d 1  *  -  o f  3  ' 

•  Transmitting  1  e  ns  -  V  9  5  ;  :m 

•  a  i  s  t  d  i  a  mo  ter  =  '3  :  -■  j  m 

•  b  r  i  :i  ;'e  s  pa  c  i ;  =  13.  b  m 

•  Collection  F"  =  5 

The  hiqh  voitu/e  to  the  phot  nmu  J  t  i  p  1  i /■  r  war- 

established  with  a  drop  lot  si  rite  of  dl  ame  l«.*  r  i  ; 

Thu  size  ranne  of  interest  was  9  to  9 1 

A  spray  of  droplets  with  prirury  diameter  of 
33  ,a:\  was  then  produced  w-inp  the  Ber:lund  I. in 

with  a  flow  rate  ot  ’.dl  cc.  mi:,  and  frequency  of 

93.dK  Hz  and  the*  dispersion  air.  For  tin-  condi¬ 
tions  selected  above,  this  size  will  produce  i 
visibility  of  3b-. 

Notice  that  for  primary  droplets  of  3  3  the 
droplets  will  have  a  diameter  of  i 1  *  \  t y) 

=  bb  and  the  triplets'  diameter  will  be  31 

/.  33  =  7b  .im . 

}  .a.  s  it-  ,.:own  for  both  Vi  s  i  hi  1  i  t  y  I  n  t  en- 

i  .  and  visibility  only  to  provide  some  compar- 

i  son . 

FijMir>.*  b  shows  the  mnis'in*m-:it.s  of  the  mono- 

d  i  s pe  r  s >■(]  s  p ray  of  5  3  ,n . 

Fi/ur.  7a  shows  the  data  of  the  spray  formed 

of  primary  droplets  u-d  doublets. 

Figure  7b  shows  results  similar  t.  >  Fi.uire  7a 

ie  i  r  obtained  witii  v  i  s  i  b  i  i  i  t  y  only.  N  » t  i  c  e  l  i .  cw 
iDurii  broad'T  the  distributions  ar*  . 

Fipnre  7  c  shows  a  spray  coot  i  i  n  i  :r..  primary 
dr«»pl«*ts,  doublets  and  triplets. 

1-  i  r  *•  7d  shows  results  sin.iiar  t  /  hl.-uire  7. 

but  ohtaineri  witn  visibility  only.  Notice  r.  hat 
the  distribution  is  broad*  t  ml  the  i  *-so  1  it  i  on  is 
not.  as  line. 

Both  the  ut-uru<v  and  r»soluti‘iu  ef  {  !u*se 
measurement"*  are  vrv  /.ood.  I  In-  t  Ivor  **t  i  .•  a  1  1  v 

,*r*-di.  ted  size*-,  are  N  i  .:m  ,  n‘-  .im,  ami  7c  ..  .  Pi*- 
c.  .r  respond  i  ;i,;  :'eanur<d  *1  i  a  n.  •  t  •  •  r  s  ire  (  [  . ,  >  )  > 

(  c.  J  to  bo  )  ,  and  I  3  ;i:;i . 

i.ar  ;ej  err  or-.-  (  ih*»ut  '•)'  )  -  n  be  e./p.-c  f  *•  1  when 

meas-irio/  tie*  s'/*-  at  a  h  ,  v  I  s  i  b  i  1  i  t  y  'Si  ). 

'  1  t  .  i 

11:1*  s  i  :iim  |  f  i : i» •  o' i *«  •  I  . '  i  r  a*  visihilitv  and 

intensify  to  :ne  i ■  i  r  *  and  v  il  »■  f  ti"  diaisr.-r  of 

Iropl.-ts  in  \  s-i  r  ty  i  >.  a*.w:  i  to  h*  i ;  .  ■  1 1  a !  ••  j-i  i 
wi  t  h  h  i  ,'h  resol  :t  i  m.  I  hi  .  :  .-i  r  n  i  *;m<*  will  <»  ■  od 

the  use  ot  Mie  j  i  t  e  . ,  a  .  i  ■.  i  *.  J  J  1  I  y  (  e  j,,,  {:;  =  ;•  »  * 

:>■•  i  r  •  r  ea  ]  j  ‘  i  ••  n  y  i  r  .  ;:*•■  o  •  -  1  will  .  *  ■  - 1  1 


erroneous  data  by  beinq  selective  duriu,.*  the  data 
validation  process.  Additional  work  is,  however, 
needed  to  establish  the  characteristics  of  the 
rejected  data. 

Since  many  of  the  erroneous  data  are  the 
result  ot  hiyh  droplet  concent  rat  i  on ,  the  percent 
>f  signals  rejected  can  he  used  to  establish  the 
limits  of  tli  technique. 

•\n  automatic  calibration  lo>»ie  used  to  estab¬ 
lish  the  ,'aie  of  the  photodetector  makes  the  lech- 
ni  quo  se  1  f'-ca  l  i  brat  i  u y . 
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Ahst  ract 

A  nonlntrusi ve  particle  sizing  technique  for 
measuring  the  size  and  velocity  of  particulates  In 
a  rocket  plume  Is  presented.  This  technique 
utilizes  two  crossing  plane  laser  sheets  to  form 
an  elongated  probe  volume  wl(h  a  known  Intensity 
profile.  This  shape  of  probe  volume  is  not  only 
less  susceptible  to  beam  wander,  but  also  has  the 
capability  of  two  dimensional  probing  of  particle 
size  and  velocity.  Experimental  results  are  pre¬ 
sented  for  particles  in  the  size  range  of  1.1  pm 
to  3.3  pm.  These  results  show  good  consistency 
with  Mie  calculations.  Actual  measurements  on 
rocket  engine  plumes  would  require  a  more  powerful 
laser  and  faster  electronic  processors. 

lot  roduct  ton 

The  mechanism  of  particle  formation  In  the 
combustor  and  plume  of  a  rocket  engine  is  not  well 
understood.'  These  particulates  can  cause  nozzle 
erosion  resulting  in  performance  degradation.  A 
tungsten-tipped  ablatively  cooled  sampling  probe, 
immersed  In  a  high-temperature  (~AOOOK)  environ¬ 
ment  with  a  large  concentration  of  high-velocity 
aluminum  oxide  particles  from  the  rocket  motor 
exhaust,  revealed  that  the  solid  particulates  of 
were  of  submicron  size  and  formed  agglomer¬ 
ates  of  up  to  a  few  microns.1  Thus,  the  ability 
to  measure  small  particles  (submicron)  moving  at 
high  speed  ('10  m/s)  In  an  optically  noisy 
environment  Is  the  key  to  rocket  engine  research 
and  performance  assessment. 

A  laser  plane  veloclmetry  technique  is  pre¬ 
sented  In  this  paper.  It  Is  based  on  laser 
Doppler  veloclmetry,  but  laser  beams  are  expanded 
into  plane  sheets  which  cross  to  form  an  elongated 
probe  volume.  Parametric  numerical  calculations 
if  the  Mie  scattering  solution  were  performed  to 
establish  the  optical  conditions  less  subject  to 
error  and  ambiguity.  An  optical  breadboard  was 
built  based  on  the  results  of  the  above  studies. 
This  breadboard  was  Interfaced  to  an  electronic 
data  processor  and  a  data  nwnagement  system. 

Known  size  polystyrene  mi crospheres  were  used 
to  calibrate  the  system.  Arbitrary  distributions 
were  formed  by  mixing  different  size  polystyrene 
particles.  Results  are  presented  for  monodls- 
perse,  blmodal  and  trlmodal  distributions.  A 
Berglund-Llu  droplet  generator  which  is  capable  of 
producing  strings  of  known  size  droplets2  (>40pm) 
was  also  used.  Smaller  size  particles  were  slau- 
lated  by  masking  the  receiver.  This  method  was 
primarily  used  to  align  the  optics  and  estimate 
the  signal  to  noise  ratio.  Studies  of  the  signal 
to  noise  ratio  as  a  function  of  particle  size  at 
various  photoimi It i pi ler  gains  showed  that  a  higher 
power  laser  would  be  required  to  extend  the 
measurable  range  to  smaller  size  particles.  For 
an  actual  rocket  engine  plume  measurement,  a 
faster  electronic  processor  would  also  be  needed. 
CopyritM  ©  America*  ImlilaK  of  Acroaaatks  lad 
Aitramtin,  Ik.,  IMS.  AN  rights  reserved. 


Background 

For  particles  significantly  larger  than  Che 
wavelength  most  of  the  light  is  scattered  Into  the 
forward  direction.  The  Fraunhofer  diffraction 
theory  adequately  describes  this  scattering  pheno¬ 
menon.3  There  are  several  particle  sizing  tech¬ 
niques  based  on  this  principle.  Swithenbank  et 
al1*  utilized  a  monolithic  photodiode  array 
detector  to  analyze  the  spatial  distribution  of 
he  scattered  light.  They  then  obtained  the  size 
distribution  parameters  form  from  Che  Fourier 
transform  of  the  signal  amplitude  distributions. 
For  particles  in  the  size  range  of  interest,  the 
scattering  function  is  a  very  complex  function  of 
the  size,  shape.  Index  of  refraction  of  the  parti¬ 
cle  and  the  scattered  angle.3  The  scattered  light 
of  large  particles  can  be  described  by  three  com¬ 
ponents;  namely,  reflection,  refraction  and  dif¬ 
fraction.  The  first  two  components  are  sensitive 
functions  of  the  above  parameters.  On  the  other 
hand,  the  diffracted  light  Intensity  Is  indepen¬ 
dent  of  the  refractive  index  and  monotonies  lly 
Increases  with  size.  Since  the  diffraction  pat¬ 
tern  of  the  scattered  light  Is  concentrated  in  the 
forward  direction.  It  leads  to  a  useful  concept 
for  particle  sizing.  The  Intensity  of  the 
scattered  light  collected  in  the  forward  direction 
will  be  Insensitive  to  the  variation  of  the 
refractive  index  and  particle  shape  and  be  a 

monotontcally  Increasing  function  of  particle 
size.  However,  a  forward  scatter  geometry  implies 
a  very  large  measurable  volume.  Moreover,  beam 
stops  are  required  to  block  the  unscattered  laser 
beams  which  can  generate  noise  to  the  receiver.  A 
reasonable  compromise  would  be  a  near  forward 
scattered  angle  of  collection  to  avoid  the  inci¬ 
dent  laser  beams,  and  reduce  the  size  of  the  probe 
volume.  It  should  be  noted  that  for  the  size 

range  of  interest  In  rocket  plumes  (0.3  pm  to 

10  pm),  Fraunhofer  diffraction  cannot  be  used. 
However,  we  found  by  solving  the  rigorous  Mie 
scattering  equations  that  many  of  the  desirable 
features  of  the  forward  scattered  light  were  also 
encountered  at  that  size  range. 

Several  optical  scattering  techniques  have 
been  proposed  for  measuring  particle  size  and 
velocity.  Very  few,  however,  would  work  In  a 

rocket  plume  environment.  The  visibility  tech¬ 
nique5,  for  Instance,  Is  only  applicable  to  parti¬ 
cles  imich  larger  than  the  wavelength  and,  as  was 
pointed  out6,  beam  degradation  will  introduce  very 
large  errors.  The  intensity  ratloing  technique7 
has  been  applied  to  small  particles  In  a  more 
limited  size  range  (.3  pm  to  2  pm). 

IMAX/flMAX8  which  use  three  laser  beams 
crossing  at  the  probe  volume  are  capable  of 
measuring  particles  in  the  size  range  of  Inter¬ 
est.  However,  it  is  expected  that  at  the  rocket 
engine  nozzle  expansion  region,  large  temperature 
gradients  could  cause  beam  wander  which  will  limit 
the  usability  of  these  techniques.  A  technique 
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based  on  Che  scattered  light  intensity  has  demon¬ 
strated  capability  of  measuring  particle  size  dis¬ 
tributions  in  the  1-30  um  range  in  a  flow  of  high 
temperature.^  We  extend  this  method  to  a  plane 
laser  Doppler  velocimetrv  type  of  measurement 
where  particles  are  illuminated  with  a  known  light 
intensity.  Particle  sizes  can  be  obtained  from 
the  scattered  light  intensity  compared  to  Mie 
calculations.  This  method  requires  a  known  size 
particle  to  calibrate  the  system. 

Description  of  Technique 


a  series  of  regions,  each  one  with  uniform, 
although  different,  Intensities.  Each  one  of 
these  regions  can  be  imaged  to  a  photodetector, 
thus  providing  a  simultaneous  measurement  of  the 
particle  size  and  velocity  at  different  parts  of 
the  spray. 

This  is  best  illustrated  with  the  mathematical 
expressions  that  define  the  above  concept.  The 
intensity  o(j  the  light  scattered  r>v  a  particle  of 

size  a  =  -  ,  index  of  refraction  n,  at  an  angle 

9  is  given  by: 


In  this  technique,  laser  sheets  are  formed 
with  cylindrical  lenses  and  crossed  to  form  a 
probe  volume  (Figure  1).  This  produces  an  elon¬ 
gated  probe  volume  which  provides  the  potential  of 
measuring  various  positions  of  the  plume  at  the 
same  time.  The  size  of  the  particles  is  obtained 
from  the  light  scattered  when  they  cross  the  probe 
volume,  and  the  velocity,  from  the  Doppler  fre¬ 
quency.  Because  the  probe  volume  is  elongated, 
the  intensity  distribution  is  a  very  "flat" 
Gaussian,  and  we  can  think  of  the  probe  volume  as 


(b) 

Fig.  1(a)  F.xpanded  laser  sheets  of  wavelength  X 
crossing  at  angle  y  to  form  an  elon¬ 
gated  probe  volume.  Particles  move 
along  x  as  shown. 

(b)  Light  intensity  at  each  portion  of  the 
probe  volume  r0,rj...is  known  from  the 
elongated  Gaussian  profile.  Particle 
size  can  be  determined  from  the  light 
scattered  from  each  portion  of  the 
probe  volume.  6  is  the  fringe  spacing. 
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where  d  is  the  particle  diameter,  V  is  its  visi¬ 
bility,  X  the  wavelength,  K  the  scattering 
coefficient,  y  the  cross-section  angle,  b(j  the 
waist  radius  in  the  x  direction,  bj  the  waist 
radius  in  the  y  direction,  and  x,v,z  the  coordi¬ 
nates.  Origin  is  defined  at  the  center  of  the 
probe  volume. 


The  scattering  coefficient  K(u,n,9)  is,  in 
general,  a  very  complex  function  when  d  is  of  the 
order  of  X.  However,  a  parametric  s*-udv,  dis¬ 
cussed  later,  showed  that  for  shallow  angles  (9  < 
S°)  the  function  is  close  to  monotonic  with  parti¬ 
cle  size  thus  .allowing  the  measurement  of  a.  The 
particle  size  is  obtained  from  the  peak  of  the 
scattered  light  after  the  ac  modulation  is  fil¬ 
tered  out.  This  function  is  normally  referred  to 
as  pedestal,  and  is  given  by 

2  2  2  2 
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The  peak,  of  the  scattered  light  for  particles 

traveling  in  the  x  direction  occurs  at  x  =  0. 

Also  the  variation  of  the  intensity  in  the  z 

direction  can  he  normally  neglected.  Therefore, 
the  peak  of  the  pedestal  at  any  location  y  is 

g  i  ve  n  by 


P  (y)  =  21  K  (a,n,9)  exp  r  .  (!) 

max  o  ,  2 1 

bl 

Since  the  laser  beanw  are  expanded  in  the  y  direc¬ 
tion  by  cylindrical  lenses,  bj  can  be  quite  large 
(typically  bj  *  20  b,j).  The  probe  volume  can  then 
be  divided  into  various  segments  of  width  Ay,  such 
that  each  one  of  these  segments  has  an  associated 
value  of  y  and  for  all  practical  purposes  the 
intensity  within  any  such  segment  is  constant. 
Therefore , 

"> 

P  (v  )  =  21  K(a,n,9)  exp  f  -  2(1  1  (O 

max  i  o  '  ,  2 


hacn  of  these  phot  ocletert  ors  wi i 1  measure 

P  (v  -  ).  and  knowing  i  and  Av  one  can  solve  lor 
max'7!*  * 

K(a,n,9).  For  known  values  of  0  and  n  one  can 
then  solve  for  the  size  parameter  a  =  - —  and  from 
it  obtain  the  diameter  of  the  particle. 
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Comput it  ions 


Scat  ter i ng 


To  obtain  the  size  parameter  a  it  is  necessary 
to  know  the  functional  relationship  of  KCajn.O)  as 
niven  hv  Ka nation  (4).  As  mentioned  earlier,  this 


iat  ions.  Both  real  ( n  j )  and  imaginary  (n->) 
of  the  refractive  indices  have  been  vari 
check  the  sensitivity  of  these  parameters, 
values  were  obtained  from  the  document  AFR 
•Sl-'jA.  Figures  2  and  3  show  the  sea 
intensity  as  i  function  of  the  particle 


•Sl-'jA.  Figures  2  and  3  show  the  sea 
intensity  as  i  function  of  the  particle 
parameter  a  for  different  values  of  n.  and 
different  scattered  angles.  Figure  4  stu 
integrated  scattered  intensity  for  a  F/d.h  c 
t i on  lens  centered  at  different  scattered  anj 

Tr  f'fMii'  1 1 fhptjo  i'j  I  nil 


DIAMETER  (vim) 


DIAMETER  (yin) 

Integrated  intensity  function  for  0=5° 
and  F/5  lens  for  polystyrene  micro¬ 
spheres,  (n  =  1.S9)  as  a  function  of 

diameter. 


Fig,  4  Integrated  intensity  function  for  a 

collection  lens  F/S.f>  and  relative  index 
of  refraction  n  =  1.81-11,1*)  x  10“’  at 
0  s  0°,  *)°,  10°,  20°. 

(dnj/dT)  is  about  10”^  1/K11  and  therefore,  we 

should  expect  small  variations  of  tij  in  the  rocket 
plume. 

Since  the  testing  and  calibration  of  the 
system  will  be  conducted  with  known  size  poly¬ 
styrene  particles,  the  scattering  function  corre¬ 
sponding  to  these  particles  was  also  evaluated. 
The  results  .are  shown  on  Figure  *)  and,  as  before, 
the  shallow  angles  offer  the  best  conditions. 

Experimental  Se  tup 

Figure  6  shows  a  schematic  diagram  of  a  system 
using  plane  laser  sheets  to  measure  the  size  and 
velocity  of  particles  in  a  flow.  Mirrors  Ml  and 
M2  steer  a  He-Ve  laser  beam  Into  a  beam  split  ter 
BS.  The  separation  of  the  two  beams  from  the  beam 


splitter  can  be  adjusted  by  changing  the  position 
of  the  input  laser  beam.  By  rotating  the  lx?  am 
splitter  the  r.ngle  between  the  two  beams  can  be 
adjusted.  Two  cylindrical  lenses  CLl  and  CL2 
expanded  one  dimension  of  the  laser  beams  into 
plane  sheets.  Another  cylindrical  lens  CL3 
focused  and  crossed  the  two  laser  sheets  at  their 
beam  waists.  The  two  laser  beams  were  terminated 
by  beam  stoppers.  Light  scattered  from  particles 
crossing  the  probe  volume  was  collected  at  *)°  by 
lenses  SLl  and  Sl,2  of  focal  lengths  K)*>  mm  and  762 
mm,  respectively*  The  collected  scattered  light 
was  focused  on  a  400  ym  di amt1  ter  pinhole  in  front 
of  a  photomultiplier  tube  (EMI  9781R)  thus 
limiting  the  field  of  view  of  the  probe  volume. 
This  insures  a  constant  laser  intensity  along  the 
slit  width.  The  f^fT  was  then  connected  to  a 
VPlpOl  signal  processor  which  was  interfaced  to  a 
computer.  The  velocity  of  the  particles  was 
measured  from  the  frequency  of  Doppler  bursts. 
The  size  of  the  particles  was  obtained  form  the 
peak  of  the  absolute  scattered  light  using  a  pulse 
height  analyzer  interfaced  to  the  frequency 
counter.  A  criterion  was  also  set  hv  the  signal 


Fig,  6  Schematic  Representation  of  Experiment. 

processor  to  reject  optical  noise  and  to  limit  the 
detection  of  the  scattered  light  to  particles  in 
the  probe  volume.  A  low  pass  filter  eliminated 
the  ac  component  of  the  signals  measured  by  the 
pulse  height  analyzer. 


Results 


Two  approaches  were  used  to  produce  known  size 
particles  to  calibrate  the  system  and  estimate  the 
S/N.  The  first  consisted  of  using  a  Berglund-Liu 
monodisperse  droplet  generator  which  capable  of 
producing  a  string  of  equal  size  dr:>.  lets.  Figure 
7  shows  a  typical  size/velocf'  ,tograra  of  this 
monodisperse  distribution.  T»  droplet  size  in 
this  case  is  41  pm,  which  ich  larger  than  the 
size  range  of  interest  *  rocket  plume  diagnos¬ 
tics.  Since,  these  d’  >plets  are  large  (40  to  60 
pm),  they  were  only  used  to  align  the  optics  and 
obtain  some  simul ±ced  conditions  by  masking  the 
receiver  and  limiting  the  collected  scattered 
light.  Thus,  particles  of  different  sizes  could 
be  simulated  and  the  S/N  established. 
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Fig.  7  Monodisperse  droplets  from  Berglund-Liu 
droplet  generator  of  diameter  of  40.8  pm. 

The  S/N  study  was  conducted  at  various  gain 
and  threshold  levels  yielding  the  conditions  where 
satisfactory  measurements  could  be  made.  Note 
that  the  detection  criterion  is  based  on  the 
ability  of  the  ac  signal  to  produce  N  fringes 
above  the  threshold  with  a  frequency  stable  to 


within  an  error  established  by  the  electronics. 
Normally,  this  ac  signal  results  from  the  fringes 
which  limit  the  probe  volume.  It  was,  however, 
found  that  for  signals  with  large  optical  noise, 
the  light  scattered  from  a  single  beam  could  pro¬ 
duce  sufficient  ac  content  that  would  fool  the 
electronics  and  be  processed.  Figure  8  shows  the 
threshold  as  a  function  of  diameter  for  different 
high  voltages  to  the  photomultiplier.  The  dia¬ 
meter  shown  is  the  abscissa  normalized  to  unity 
and  represents  the  relative  position  within  the 
size  histogram.  This  result  indicates  that  an 
intense  light  source  and  a  low  PMT  voltage  will  hj 
required  to  obtain  a  large  dynamic  size  range 
given  a  specified  threshold  level. 


Fig.  8  Noise  level  as  a  function  of  particle 
size  for  different  PMT  voltages. 


The  second  approach  consisted  of  blowing 
polystyrene  particles  of  1.1  pm  to  3.3  pm  in 
diameter.  These  particles  were  first  diluted  in 
distilled  water.  A  nebulizer  powered  by  an  air 
compressor  was  then  used  to  inject  the  water  mist 
together  with  the  polystyrene  particles  into  a 
heating  chamber.  There  the  water  mist  was  evapo¬ 
rated  and  the  solid  polystyrene  particles  were 
blown  (by  the  same  air  compressor)  into  the  prob-' 
volume.  This  allowed  uie  measurement  of  known 
size  particles  in  the  size  range  of  interest. 
Different  size  particles  could  be  diluted  in  water 
thus  producing  more  complex  distributions  such  as 
hi  modal  and  trimodal. 

Figure  9  shows  typical  histograms  (size  and 
velocity)  of  the  1.7  pm  particles.  Figure  10 
shows  the  results  corresponding  to  a  bimodal 
distribution  of  l.l  pm  and  1.7  pm,  and  Figure  11 
shows  a  trimodal  distribution  of  1.7  pm,  2.7  pm, 
and  3.3  pm.  It  should  be  noted  that  the  high 
voltages  to  the  PMT  in  Figures  10  and  11  are 
different  and,  therefore,  the  location  of  the  1.7 
pm  is  shifted.  We  had  to  do  this  because,  in 
order  to  measure  the  l.l  pm,  the  high  voltage  had 
to  be  Increased.  This  high  voltage  was,  however, 
not  compatible  with  the  3.3  pm  particles  since, 
for  the  limited  laser  power  used  in  these  exper¬ 
iments  (23  nW  of  red  Light),  the  3.3  pm  signal 
would  have  been  much  broader  due  to  excess  noise. 

The  theoretical  relative  sizes  obtained  from, 
the  Mie  codes  are  shown  on  the  histograms.  For 
the  most  part,  the  peaks  are  within  10?.  of  the 
predicted  ones  which  is  consistent  with  sow  of 


V 


the  ambiguities  experienced  by  the  scattering 
function.  it  can  be  observed  in  the  histograms 
that  the  size  dis t r i but i on  of  the  larger  particles 
Is  broader.  This  broadening  would  be  considerably 
reduced  with  a  more  powerful  laser  or  more 
advanced  electronic  processing.  Nevertheless,  the 
results  are  quite  remarkable  for  a  noni nt  nisi ve 
technique  making  measurements  in  such  complex  size 
range.  These  results  confirm  that  the  technique 
is  extremely  promising  to  perform  measurements  in 
a  real  rocket  plume,  given  the  aval' Lability  of  a 
high  power  laser  and  a  faster  electronic 
processor. 
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Histograms  of  monodisperse  polystyrene 
microspheres  of  1.7  pm  in  diameter. 
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Fig.  10  Histograms  of  bimodaL  polystyrene  micro- 
spheres  of  l.l  pm  and  1.74  pm  in  diameter 
(arrows  indicate  the  relative  theoretical 
diameters ). 
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